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Forward-backward correlation strength (6) as a function of pesudorapidity intervals for 
experimental data from p+p non-singly diffractive collisions are compared to PYTHIA 
and PHOJET model calculations. The correlations are discussed as a function of ra- 
pidity window (Ar;) symmetric about the central rapidity as well as rapidity window 
separated by a gap (r) gap ) between forward and backward regions. While the correla- 
tions are observed to be independent of Ar], it is found to decrease with increase in r) gap . 
This reflects the role of short range correlations and justifies the use of r) gav to obtain 
the accurate information about the physics of interest, the long range correlations. The 
correlation strength from PYTHIA are in agreement with the available experimental 
data while those from the PHOJET give higher values. For p+p collisions at \fs = 7, 
10 and 14 TeV, the correlation strength from PHOJET are lower compared to those 
from PYTHIA, this is in contrast to the observations at lower energies. The experimen- 
tal b value shows a linear dependence on In yTs with the maximum value of unity being 
reached at ^/s = 16 TeV, beyond the top LHC energy. However calculations from the 
PYTHIA and PHOJET models indicate a deviation from linear dependence on lny^s 
and saturation in the b values being reached beyond ^/s = 1.8 TeV. Such a saturation in 
correlation values could have interesting physical interpretations related to clan struc- 
tures in particle production. Strong forward-backward correlations are associated with 
cluster production in the collisions. The average number of charged particles to which 
the clusters fragments, called the cluster size, are found to also increase linearly with 
In y/s for both data and the models studied. The rate of increase in cluster size vs. lny^s 
from models studied are larger compared to those from the data and higher for PHOJET 
compared to PYTHIA. Our study indicates that the forward-backward measurements 
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will provide a clear distinguishing observable for the models studied at LHC energies. 

Keywords: Forward-Backward correlations, Long range correlations, cluster formation, 
p-|-p collisions at LHC 

PACS numbers: 1.200, 1.300 

1. Introduction 

Understanding the mechanism of particle production in p+p collisions is one of 
the first goals of the experiments at the Large Hadron Collider (LHC). Several 
new results from the LHC experiments in terms of particle multiplicities at y/s — 
0.9, 2.36, 2.76 and 7 TeV are being compared to various models of particle produc- 
tkfflJ ^ I 4 | 5|6 | 7 | The most popular model s for co mparison to data from p+p collisions 
being PYTHIA I 8 ! 9 ' 10 ! 11 ! and PHO JET U M2i T he PHO JET model combines the 
ideas based on a dual parton model (DPM) I 14 | 15 | on soft process of particle produc- 
tion and uses lowest-order perturbative QCD for hard process. Regge phenomenol- 
ogy is used to parameterize the total, elastic and inelastic cross-sections. The initial 
and final state parton shower are generated in leading log-approximation. PYTHIA 
on the other hand uses string fragmentation as a process of hadronization and tends 
to use the perturbative parton-parton scattering for low to high pt particle pro- 
duction. Initial results have shown that both models do not ha ve perfect a greement 
with multiplicity measurements at LHC energies studied so far HEHEEEI, However 
it may be mentioned that several of these models in turn are used to obtain various 
correction factors for the experimental measurements. In this paper we suggest that 
the correlations between the particles produced in forward and backward rapidities 
can be used to discriminate between the various models of particle production, more 
reliably at the LHC energies. 

The forward-backward correlations previously observed had several physical in- 
terpretations. The correlations over small range in rapidity are believed to be dom- 
inated by short-range correlations as due to resonance decays, those occurring over 
large rapidity range could be interpreted to be due to multiple parton interactions. 
In early 1985, a statistical interpretation of the forward-backward correlations ob- 
served in ISR energies was provided ^1 The interpretation was based on clusters 
being produced in these hadronic collisions according to a negative binomial distri- 
bution and the final hadrons are a result of the decay of thes e cl usters. An extended 
version of such a statistical scenario can be found in Ref. ^[ A more dynamical 
interpretation was provided based on the DPM as in Ref. The experimental data 
was also interpreted in terms of a model based on a unitarized model which included 
soft and semi-hard components (minijets). In such a model the average number of 
particles produced was proportional to the effective number of inelastic collisions. 
The increase of forward-backward correlations with yfs was fully generated by the 
superp osit ion of the different impact parameter contributions to the inelastic cross 
section ^1. The extension of this interesting idea is to consider a weighted superpo- 
sition of two classes of events in hadronic collisions: soft and semi-hard processes. 
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The behavior of the semi-hard component on observed forward-backward correla- 
tions can lead to interesting interpretation in terms of the clan structure of particle 
production as discussed in Refs ^1. The authors of Refs ^1 envision a possibility of 
formation of new species of clans at LHC and a possible phase transition in clan 
production mechanism. 

Recently the forward-backward correlations have been related to the simplest 
form of partonic interaction that exhibits back-to-back correlation. Such a formula- 
tion is found to give a fairly good description of data from STAR Collaboration in 
p+p collisions at ^fs = 200 GeV Studies of forward-backward correlations pro- 
vide baseline measurements to look for long range rapidity correlation in heavy-ion 
collisions. STAR experiment has recently made this relative (comparison between 
p+p and different collision centrality Au+Au collisions at the center of mass energies 
of 200 GeV) measurements to claim the existence of a large long range correlation 
in central Au+Au collisions ^21. Although the experimental data cannot differen- 
tiate whethe r the actual underlying mechanis m is d ue to those in a dual parton 
DPM or a Color Glass Condensate both however require that the 
long range correlations are produced by multiple parton-parton interactions. Fur- 
ther it is ar gued that the clustering of color sources could lead to forward-backward 
correlations^^. Hence the measurement of the long range forward backward corre- 
lations in the multiplicity of produced particles in high energy collisions will give 
us insight into the space-time dynamics of the collision. 

In this paper, we first review the existing experimental data on forward-backward 
correlations, discuss the model calculations and methods used to extract these cor- 
relations. This is followed by comparison of the experimental data on forward- 
backward correlations with those from PYTHIA and PHO JET models, expectations 
at top LHC energies and interpretation of the data in terms of cluster production in 
high energy hadronic collisions. We find that the forward-backward correlations in 
p+p collisions at y/s = 7, 10, 14 TeV will help differentiate various models and hence 
the underlying particle production mechanism. Saturation of correlation strength 
at LHC energies could indicate new physics. 

2. Experimental data, Model Simulation and Correlations 

The experimental charged particle data reviewed in this paper are from the E735 
and UA5 collaboration. These data are for p+p collisions and corresponds to non- 
singly diffractive events. The E735 data ^corresponds to = 300, 546, 1000 and 
1800 GeV while those from UA5 27 corresponds to y/s = 200, 546 and 900 GeV. 
The pesudorapidity acceptance range for the E735 experimental data is | i] \ < 3.25 
and those for UA5 is | r\ \ < 4.0. 

The model results presented in this paper are from PYTHIA (Version 6.4) and 
PHO JET (Version 1.1, uses jetset74 from PYTHIA) with default settings. The 
event type selected are non-singly diffractive as for the existing measurements. For 
comparison to existing data the simulations are done for p+p collisions with the 



March 3, 2013 7:18 WSPC/INSTRUCTION FILE fb-ijmpa'revised 



4 A. K. Dash, D. P. Mahapatra and B. Mohanty 



! r< — 


^gap 


— J*-! 1 




^— ~ ! ! 




Ar| 





-T)2 -r\l r\ = rilri2 r\ 



Fig. 1. Definition of Ar? and ?7gap used in forward- backward correlation studies. 



experimental acceptances included. For predictions at LHC energies of y/s = 10 
and 14 TeV the model calculations are done for p+p collisions. A realistic transverse 
momentum cut of above 100 MeV for charged particles are applied. 
The forward-backward correlations can be obtained by two methods. 

• If TVf, is the multiplicity in backward hemisphere and Nf is the correspond- 
ing multiplicity in forward hemisphere for the same event, then the corre- 
lation co-efficient can be obtained by plotting (Nt)(Nf) vs. Nf. Where (..) 
denotes event average. The resultant distribution can be fitted to a linear 
function as (Nb)(Nf) = a + b Nf, to obtained the correlation strength, b. 

• The correlation coefficient can also be defined as b = ^ J fl^L ^/M Nb ' > ■ 
Most of the results presented in this paper uses this method for calculating 
b. For the model results presented we have explicitly checked and found 
there is good agreement between the two methods. 

Before we proceed towards discussion of the experimental data and results from 
model calculations, we discuss few relevant experimental aspects. The forward- 
backward correlation results for a given y/s are usually presented as a function 
Ar] and r/gap- Fig. Q] shows how these quantities are defined. The A77 corresponds 
to calculation of correlation coefficient in a symmetric window about the central 
rapidity. Increasing Arj window includes the contributions from smaller Arj win- 
dow. The results as a function of 77 ga p corresponds to correlation co-efficient being 
calculated in some fixed Ar/ window separated by a gap in rapidity between forward 
and backward regions by an amount r/gap ■ 

Most of the existing forward-backward correlation results at high energies are 
from p+p collisions. The current work presents the predictions of the correlation 
co-efficient from PYHTIA and PHOJET models in p+p collisions at LHC energies 
and compares them to extrapolations of results from p+p collisions. So it is essen- 
tial to check using models if a difference is expected in correlation values between 
p+p and p+p collisions at a given y/s. Note however that the STAR Collaboration 
has measured the strength of charged particle forward-backward multiplicity cor- 
relation, b in p+p collision at y/s — 200 GeV® . It is about 3-4 times smaller 
than the one measured by UA5 in p+p collision at the same energy apparently ^3. 
The eta gap in STAR is 0.2 where as in UA5 it is 0.5. Also both the experimental 
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Fig. 2. Forward-backward correlation co-efficient from PYTHIA at s/s = 900 GeV as a function 
of Ar; and r] gap . The results are from non-singly diffractive events and presented for both p+p 
and p+p collisions. The effect of weak decay, lower transverse momentum cut off of 100 MeV/c, 
typical charged particle detector efficiency, different model parameter are studied. For comparison 
the results from p+p collisions in UA5 experiment are also shown. The bottom panels of the figure 
shows the ratio of the correlations from simulations with various conditions to those measured by 
UA5 experiment. 



acceptances are different. STAR has | 77 | < 1.0 where as UA5 has | rj \ < 4. 
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Fig.[5]shows the correlation co-efficient (b) calculated using non-singly diffractive 
events from PYTHIA model at </s = 900 GeV as a function of Aij for p+p and p+p 
collisions. No appreciable difference is observed between the two colliding systems for 
correlations studied as a function of Ar] and r] gap . Also shown for comparison are the 
charged particle correlation co-efficient results from UA5 experiment. Usually due to 
experimental limitations the analysis is carried out with a lower pt cut off. As seen 
in Fig. [2] a typical experimental cut-off of 100 MeV on pt does not seem to change 
the correlation values as a function of Arj and r] gap . We have also investigated the 
effect of weak decay particle contributions to this analysis. The difference between 
weak decay on and weak decay off for p+p collisions is noticeable for smaller Ar] 
and the difference vanishes as we go to larger Ar]. This is along expected lines, 
as decay effect will introduce short-range correlations. The effect of finite charged 
particle detection efficiency has been studied by varying it between 86% to 94%, no 
appericiable change is observed. The effect of a different tuned version of PYTHIA 
(D6T) has also been investigated. This also does not seem to affect the observed 
correlations by a large amount. The bottom panel of the figure shows the ratio of 
the correlation from simulation with various effects to the correlation measured by 
UA5 experiment. One observes that for the study with respect to Arj the variations 
lie within 20%. For the study with respect to the rj gap , the ratio is close to unity 
upto a value of r] gap 4 units, then the simulation starts to diverge away from the 
measurements. 

3. Correlations in Ar] and ?7 gap 

Fig [3] shows the correlation strength (6) as a function of At? for non-singly diffrac- 
tive events in p+p collisions from E735 ^1 anc j UA5 experiments at various y/s 
compared to corresponding results from PYTHIA and PHOJET model calculations. 
The following observations can be made: (a) The correlation strength both in ex- 
perimental data and simulations are observed to be almost independent of Ary, (b) 
correlations seem to increase with y/s, and (c) PYTHIA model calculations agrees 
well with the experimental data, while PHOJET tends to give higher correlations 
for all the measured y/s. This can be seen from the ratio of the correlation value 
from model to data in the bottom panel of the Fig. [3] As the larger Ar] includes 
contribution from smaller At] intervals, such a data may not provide accurate infor- 
mation about long range correlations. These correlations could include significant 
short-range correlations. 

Fig |4] shows the correlation strength (b) as a function of ry ga p. The value of b 
decreases with increase in 77 gap indicating the diminishing contribution from short- 
range correlations. With a 77 gap of around 2 units still significant correlation at the 
level of 30 - 40% are observed. Both the models reproduce the decreasing trend 
of the correlations with increasing %j, p . PYTHIA model calculations has better 
agreement with the measurements, while PHOJET model in general over estimates 
the correlation strength at most of the y/s, except for y/s = 1.8 TeV. This can be 



March 3, 2013 7:18 WSPC/INSTRUCTION FILE fb-ijmpa'revised 



Expectation of forward-backward rapidity correlations inp+p collisions at the LHC energies 7 



1 

pp NSD 

E735 PYTHIA 

* 1B00 GeV - ■ 1800 GeV 

A 1000 GeV 1000 GeV 

■ 546 GeV --■ 546 GeV 

• 300 GeV — 300 GeV 

I | I I I I I I I I 


i | i | i | i | i | i | i | i | i 

— ♦ — * 

UA5 PYTHIA 

A 900 GeV 900 GeV 

■ 546 GeV ■ ■ ■ 546 GeV 
- • 200 GeV — 200 GeV 


# ♦ * • 

E735 PHOJET 

- * 1800 GeV - ■ ■ 1800 GeV 

A 1000 GeV 1000 GeV 

■ 546 GeV 546 GeV 

- • 300 GeV — 300 GeV 

I 


I I I I I I I I I 

• * * * • 

UA5 PHOJET 

A 900 GeV 900 GeV 

■ 546 GeV ---546 GeV 
• 200 GeV —200 GeV 

,1,1,1,1,1,1,1,1, 



U ' u 0.5 1 1.5 2 2.5 3 3.5 4 4} 0.5 1 1 .5 2 2.5 3 3.5 4 4.5 



At] At] 



1.4 
1.2 

q 

(01.0 
0.8 
0.6 
1.4 
1.2 

g 

'ra 1.0 
tr 

0.8 
0.6 



I 1 1 1 1 1 1 1 I I I 1 I I I I 

' ppNSD 


I i I i I i I i I i I i I i I 
PYTHIA /UA5 DATA " 
900 GeV 
"* "■546 GeV 
-200 GeV 


PYTHIA /E735 DATA 
- ■■■1800 GeV 

1 000 GeV 

— 546 GeV 
" - 300 GeV . . " 


I I I I I I I 


| I | I | I | I I I I 




PHOJET / E735 DATA 
- ■1800 GeV 

- 1 000 GeV 

■■■546 GeV 

- — 300 GeV 
I 


PHOJET /UA5 DATA - 

900 GeV 

- —546 GeV 
—200 GeV 





0.5 1 1.5 2 2.5 3 3.5 4 4J 0.5 1 1.5 2 2.5 3 3.5 4 4.5 



Fig. 3. Top panel: Forward-backward correlation strength (ft) as a function of An for p+p non- 
singly diffractive collisions at various y/s from E735 and UA5 experiments. The measurements 
are compared to PYTHIA and PHOJET model calculations. Bottom panel: The ratio of model 
calculations to experimental data. 



seen from the ratio of correlations from models to that from data as shown in the 
bottom panel of Fig. S) The agreement with the highest energy data by both the 
models poses the question whether at LHC energies the current measurements will 
have the distinguishing power. To investigate this possibility, we now discuss the 
predictions at three LHC energies of \fs — 7, 10 and 14 TeV from PYTHIA and 
PHOJET models. 

Fig. [5] and [6] shows the prediction of forward-backward correlation strength (b) 
from PYTHIA and PHOJET models as a function of Arj and 77 gap for p+p collisions 
at yfs = 7, 10 and 14 TeV. Also shown for understanding are the energy dependence 
trend of the existing correlation measurements from p+p collisions. In general both 
models predict a higher but constant value of 6 as a function of A77 at y/s = 7, 10 
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Fig. 4. Top panel: Forward-backward correlation strength (6) as a function of r? gap for p+p non- 
singly diffractive collisions at various ^/s from E735 and UA5 experiments. The measurements 
are compared to PYTHIA and PHOJET model calculations. Bottom panel: The ratio of model 
calculations to experimental data. 



and 14 TcV p+p collisions, similar to the trend seen at lower energies p+p collisions. 
Although PYTHIA tends to predict a slightly increasing trend at larger Arj. Both 
models also predict a decreasing trend of b as a function of 77 gap at \fs = 7, 10 
and 14 TeV as was observed for lower energies. The decrease in b with respect to 
% ap horn PYTHIA is slower compared to that from PHOJET at 7, 10 and 14 TeV. 
Both models clearly suggest that the correlation strength should saturate at higher 
energies, with PHOJET indicating that saturation could occur as early as \/s = 1.8 
TeV. 

The energy dependence of the forward-backward correlation can be seen in the 
Fig. El All available experimental data from ISR UA5 ^3 and E735 ^ Col- 
laborations shows a linear increase in the correlation value with beam energy. The 
best description of the data is obtained as a + /31n-y/s, where a = -0.18 ± 0.02 and 
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Fig. 5. Expected forward-backward correlation strength (b) estimated from PYTHIA and PHO- 
JET models as a function of An for p+p collisions at LHC energies of \fs = 7, 10 and 14 TeV. 
Also shown for understanding are the energy dependence of the correlations from the existing 
measurements in p+p collisions. 



(3 = 0.122 ± 0.005 are the fit parameters. For the ^fs > 200 GeV the correlation 
values from PYTHIA model are in good agreement with the measurements, but 
those from PHOJET model are higher. However interestingly, at the LHC ener- 
gies the estimates for b are lower from PYTHIA compared to PHOJET. If the a 
+ Pln^/s dependence of b on y/s holds then the maximum correlation value of 1 
will be attained for y/s <~ 16 TeV, beyond LHC energies. Beyond which the value 
is expected to saturate. But the model estimates show that the values of b could 
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Fig. 6. Expected forward-backward correlation strength (b) estimated from PYTHIA and PHO- 
JET models as a function of r; gap for p+p collisions at LHC energies of y/s = 7, 10 and 14 TeV. 
Also shown for understanding are the energy dependence of the correlations for the existing mea- 
surements in p+p collisions. 



saturate or drop earlier around y/s = 1.8 TeV. The PHOJET models predicting a 
saturation/drop in correlation values starting at a somewhat smaller yfs compared 
to that predicted from PYTHIA. New data at LHC energies will provide a clear 
picture. 

Saturation or drop in the energy dependence of the co rrel ation values could have 
interesting physical consequences as is discussed in Refs In this geometric pic- 
ture, the energy dependence of b is understood based on the superposition of two 
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Fig. 7. Correlation strength, b as a function of lny^ for experimental data from ISR, UA5 and 
E735 Collaborations. Also shown are comparisons from PYTHIA and PHOJET model calculations, 
including those expected at y/s = 7, 10 and 14 TeV. Solid line is a fit to the experimental data 
points, extrapolation of which indicates the correlation strength will reach a maximum value of 
unity beyond LHC energies. 



components: soft and semi-hard in p+p(p) collisions'^. This picture at higher ener- 
gies is supported by the need of two negative binomial distributions (with diffe rent 
parameters) to explain the multiplicity distribution of produced charged hadrons^l 
Further in such an approach the particle production is thought to be from a certain 
number of independent initial sources decaying to final products. The produced 
particles from each source are expected to all stay in the same hemisphere. However 
certain leakage of particles to the other hemisphere is allowed and is controlled by 
a leakage parameter in the model. Such a model predicts a saturation of correlation 

inni 

co-efficient at LHC energies -^4 The final value of b depends on whether the leak- 
age parameter increases with y/s (higher b) , constant with y/s or decreases with y/s 
(lower 6). The prediction from this model agrees well with the results from PHOJET 
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model and is lower than PYTHIA expectations. 

4. Cluster Production 

One way to interpret the forward-backward correlations is through the correlations 
between particles originating from clusters. The average number of particles origi- 
nating from clusters is called the cluster size and can be obtained by two methods. 

• Method 1: By defining an asymmetry parameter, Z = Nf — Nb, then 
(Z 2 ) = rN c h- Where r is the cluster size and N c h is total number of 
charged particles. 

• Method 2: Directly from the measurement of the correlation coefficient 
as, b = ^"^/fc+i+r ' wnere (Nch) and k are the parameters of a negative 
binomial fit to the multiplicity distributions ^Sl. 

Figure [8] shows the typical relation between (Z 2 ) and N c h for p+p collisions 
at y/s — 900 GeV from PYTHIA. The distribution is fitted by a linear function 
(dashed and solid lines) in a fixed interval of N c h to extract the cluster size,r. We 
followed this procedure to obtain cluster size from both the models PYTHIA and 
PHOJET for various v^- 

Figure M shows the ^/s dependence of the cluster size from E735^S1, UA5^3and 
ISR^[ rp ne resu it s have been compared to model calculations from PYTHIA and 
PHOJET, as well as from a theoretical calculations by Lim et al -^l. The PYTHIA 
gives similar cluster size as seen by experiments upto 1.8 TeV and PHOJET give 
similar cluster sizes as seen by experiments up to ^/s = 546 GeV. Beyond this 
energy cluster size from PHOJET is higher. The cluster size dependence on for 
the measured experimental data can be parameterized as a' + /3'hiy/s. The values 
a' and f3' for the experimental data, PYTHIA, PHOJET and from model by Lim 
et al., are given in the Table. The expectations for top LHC energies from the two 
models are also shown. PHOJET expects the average number of particles from a 
cluster to be around 5-6, while PYTHIA gives a much lower value of around 3.1. 
Thereby providing a clear observable, when compared to LHC data, will distinguish 
between the underlying mechanism of particle production in p+p collisions. 

5. Summary 

We have reviewed the existing data on forward-backward correlations in p+p(p) 
collisions. Compared these experimental measurements to model calculations from 
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Fig. 8. Variance of the asymmetry measurement ((Z 2 )) as a function of charged particle multi- 
plicity for p+p collisions at y/s = 900 GeV from PYHTIA model. The solid and dashed lines are 
a linear fit to extract the cluster sizes in various ranges. 



PYTHIA and PHOJET. It is found that the correlation strength from PYTHIA is in 
agreement with the existing measurements, while those from PHOJET give higher 
correlations. However for top LHC energies of y/~s = 7, 10, 14 TeV, the correlation 
strength from PHOJET is lower compared to PYTHIA, suggesting a transition 
at an intermediate energy accessible at LHC. The measured correlation strength 
are found to increase linearly with In y^s and extrapolation suggests it will reach 
unity (maximum value) around -y/s = 20 TeV, beyond the beam energy reach at 
LHC. However model calculations suggest the correlation values tends to saturate 
starting at ^fs — 2 TeV. If such a saturation is observed at LHC it could mean 
interesting physical consequences related to clan structures in particle production. 
We have also reviewed the existing results on a common interpretation of forward- 
backward correlations in terms of cluster production. The cluster sizes are found to 



March 3, 2013 7:18 WSPC/INSTRUCTION FILE fb-ijmpa'revised 



14 A. K. Dash, D. P. Mahapatra and B. Mohanty 



10 



8 

CD 
N 

CD 

H — > 

c/) 



O 





pp NSD 




■ 


E735 (Method 1) 


o 


• 


E735 (Method 2) 


□ 


▲ 


UA5 (Method 2) 


A 


* 


ISR (Method 2) 






oc'+p'InVs 





NSD (Method 1) 
PYTHIA pp 
PHOJET pp 
PYTHIA pp 
PHOJET pp 
Lim □ 
□ □ 




L 



10 2 , 10 3 

Vs (GeV) 



10 4 



Fig. 9. Cluster size as a function of ^/s as measured in E725, UA5 and ISR experiments. The 
method 1 and 2 are discussed in the text. Also shown are the model comparisons with the data at 
various beam energies in p+p collisions and the expectations at top LHC energies from PYTHIA 
and PHOJET. 



increase with increase in beam energy. Similar to the correlation strength, the cluster 
size from PYTHIA compares well with the existing experimental data. For higher 
energies (> 546 GeV) PHOJET gives a higher cluster size compared to PYTHIA. 
The study of cluster size from forward-backward correlations can be a very good 
discriminator for the particle production models in p+p collisions. 
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